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Abstract 

New Zealand's endemic Stewart Island Shag (Leucocarbo chalconotus) comprises two regional groups (Otago and Foveaux 
Strait) that show consistent differentiation in relative frequencies of pied versus dark-bronze morphotypes, the extent of 
facial carunculation, body size and breeding time. We used modern and ancient DNA (mitochondrial DNA control region 
one), and morphometric approaches to investigate the phylogeography and taxonomy of L. chalconotus and its closely 
related sister species, the endemic Chatham Island Shag (I. onslowi). Our analysis shows Leucocarbo shags in southern New 
Zealand comprise two well-supported clades, each containing both pied and dark-bronze morphs. However, the combined 
monophyly of these populations is not supported, with the L chalconotus Otago lineage sister to L. onslowi. Morphometric 
analysis indicates that Leucocarbo shags from Otago are larger on average than those from Foveaux Strait. Principal co- 
ordinate analysis of morphometric data showed substantial morphological differentiation between the Otago and Foveaux 
Strait clades, and L. onslowi. The phylogeographic partitioning detected within L. chalconotus is marked, and such strong 
structure is rare for phalacrocoracid species. Our phylogenetic results, together with consistent differences in relative 
proportions of plumage morphs and facial carunculation, and concordant differentiation in body size and breeding time, 
suggest several alternative evolutionary hypotheses that require further investigation to determine the level of taxonomic 
distinctiveness that best represents the L. chalconotus Otago and Foveaux Strait clades. 

Citation: Rawlence NJ, Till CE, Scofield RP, Tennyson AJD, Collins CJ, et al. (2014) Strong Phylogeographic Structure in a Sedentary Seabird, the Stewart Island 
Shag {Leucocarbo chalconotus). PLoS ONE 9(3): e90769. doi:10.1371/joumal.pone.0090769 

Editor: William J. Etges, University of Arkansas, United States of America 

Received November 30, 2013; Accepted February 4, 2014; Published March 10, 2014 

Copyright: © 2014 Rawlence et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: Funding was provided by the Royal Society of New Zealand Marsden Fund, the Department of Zoology (University of Otago) and the Allan Wilson 
Centre for Molecular Ecology and Evolution. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the 
manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* E-mail: nic.rawlence@otago.ac.nz (NJR); martyn. kennedy@otago.ac.nz (MK) 



Introduction 

The blue-eyed shags (Leucocarbo spp.) are a species-rich seabird 
clade exhibiting a circumpolar Southern-Hemisphere distribution. 
Numerous locally endemic taxa are associated with particular 
islands or island groups [1]. Within the New Zealand region there 
are six species currendy recognised: New Zealand King Shag (L. 
carunculatus), Stewart Island Shag (L. chalconotus), Chatham Island 
Shag (L. onslowi), Auckland Island Shag (L. colensoi), Campbell 
Island Shag (L. campbelli), and Bounty Island Shag (L. ranfurlyi) [2] . 

The Stewart Island Shag L. chalconotus is endemic to southern 
New Zealand, with its current range extending from the Stewart 
Island region (Foveaux Strait) north to the Otago coast (Fig. 1). 
This species is unique within its genus in that it shows substantial 
intraspecific variation in plumage colouration, with distinct pied 
and dark-bronze morphotypes (Fig. 1). Both colour morphs are 
widespread throughout the species' range, with no evidence for 
assortative mating associated with plumage phenotype [3]. 



Nevertheless, morphological and behavioural analyses by Lalas 
[4] , and Lalas and Perriman [5] , have detected consistent regional 
differentiation between Otago and Foveaux Strait populations of 
L. chalconotus (Fig. 1), primarily associated with relative frequencies 
of pied versus dark-bronze morphotypes, the extent of facial 
carunculation, body size and breeding time. They found that, (1) 
the Otago group comprises 20-30% pied morphs versus 50-60% 
in the Foveaux Strait group (Fig. 1); (2) Otago birds have 
approximately equal frequencies of small facial caruncles versus 
scattered papillae, whereas Foveaux Strait birds always have 
scattered papillae (Fig. 1); (3) morphometric analyses of museum 
specimens revealed birds from Otago were 7-14% larger than 
those from Foveaux Strait; and (4) breeding was initiated sooner in 
Otago (May-September versus September onwards in Foveaux 
Strait) [4-5]. 

The current Otago breeding populations of L. chalconotus range 
from Maukiekie Island south to Kinakina Island, whereas the 
southern breeding populations are restricted to islands in Foveaux 
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Figure 1 . Distributional and morphological data for Chatham Island {Leucocarbo onslowi) and Stewart Island (L. chalconotus) shags. A: 

Map of New Zealand showing the location of the Chatham Islands and Otago/Foveaux Strait study sites; B: Distribution of L. onslowi (green) breeding 
colonies and roosting sites. Unfilled circles represent existing but un-sampled colonies and roosting sites. Filled green stars represent Holocene fossil 
bones at sites without colonies or roosts nearby. Leucocarbo onslowi exhibits pied plumage only (white pie graph) with pronounced bright orange 
caruncles in breeding plumage (orange pie graph). C: Distribution of L chalconotus from Otago (blue) and Foveaux Strait (red) breeding colonies and 
roosting sites. Filled circles represent samples from colonies and roosting sites, while the stars represent beach wrecks found at sites without colonies 
or roosts nearby. Unfilled circles represent existing but un-sampled colonies and roosting sites. Shags from the Otago populations have 20-30% pied 
morphs (pie graphs; black: dark-bronze; white: pied) and 50%:50% small bright orange caruncles: dark to dull orange scattered papillae in breeding 
plumage (pie graph; yellow: small bright orange caruncles; grey: dark to dull orange scattered papillae) compared to 50-60% pied and dark to dull 
orange scattered papillae in breeding plumage in shags from the Foveaux Strait populations [4-5], [1 3]. In C, multiple breeding colonies and roosting 
sites are represented at the following locations (from north to south): Warrington and Karitane; Otago Peninsula including Long Beach, 
Aramoana, Otago Harbour, Taiaroa Head, Boulder Beach, Papanui Beach, Allans Beach, Wharekakahu Island and Gull Rocks; Seal Rocks and Ruapuke 
Island; Easy Harbour and Shag Rock. 
doi:1 0.1 371 /journal.pone.0090769.g001 



Strait (Fig. 1). The current non-breeding range for the Otago 
group extends from Lake Wainono south to The Sisters in the 
Catlins (Fig. 1) [5], with rare vagrants observed as far north as 
Lake Ellesmere [6]. Historically, the Otago populations were 
centered near Otago Peninsula [7-9]. During the 1980s, Otago 
breeding colonies were located at Maukiekie Island and several 
sites around Otago Peninsula with a roost site at Nugget Point [4] 
(Fig. 1). Recent decades, however, have seen substantial northward 
(e.g. roost site at Oamaru) and southward (e.g. breeding site at 



Kinakina Island, roost site at the The Sisters in the Catlins) 
expansions of this regional population [5-6] (Fig. 1). This 
expansion results from an overall population increase and is 
apparently driven by L. chalconotus' ability to readily abandon and 
establish new breeding colonies, and large numbers at some 
breeding colonies creating overflow into adjacent regions [5]. 

There is no obvious contemporary physical barrier between the 
two regional groupings of L. chalconotus populations, and it has 
been suggested that the species may have been more widespread in 
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the past. Indeed, bones attributable to Leucocarbo have been found 
in Holocene fossil sand dune deposits and early Maori archaeo- 
logical sites (ca. 1280-1450 AD) throughout the current range of 
L. chalconotus and further north along the eastern South Island, 
suggesting a former widespread distribution prior to Polynesian 
settlement of New Zealand ca. 1280 AD [10-12]. 

The taxonomic history of L. chalconotus has been complicated 
because many of its external morphological character states are 
shared with other blue-eyed shag taxa. In breeding plumage the 
combination of caruncle size and colour is diagnostic with L. 
carunculatus having a large yellow caruncle and L. onslowi having a 
pronounced bright orange caruncle [4], [13]. In L. chalconotus from 
Otago the small caruncles are usually bright orange, whereas 
when there are only scattered papillae, these are dark to dull 
orange. Likewise, the scattered papillae in L. chalconotus from 
Foveaux Strait are dark to dull orange [4], [13]. Caruncle size 
decreases and its colour fades in non-breeding plumage, so there is 
currently no reliable method for discriminating between L. 
carunculatus, L. onslowi, and the pied morphs of L. chalconotus in 
the wild [14]. L. chalconotus and L. onslowi have previously been 
regarded as either subspecies of L. carunculatus [4], [15-16] or as 
full species [1—2], [17-19]. Siegel-Causey [18] proposed osteolog- 
ical character state changes that united his New Zealand blue-eyed 
shags, but was unable to distinguish between L. chalconotus and L. 
onslowi, due to small sample size and a lack of comparative material 
(cf. Worthy [20] who found 14 osteological character state 
differences between these two taxa). Siegel-Causey [18] tentatively 
retained L. chalconotus as a distinct species, sister to a group 
containing L. onslowi, L. ranfurlyi and L. colensoi. Worthy [10] 
criticised Siegel-Causey's [18] argument for maintaining separate 
species status for L. chalconotus when L. carunculatus was not included 
in the study. 

Based on the regional differentiation detected by Lalas [4] and 
the absence of apparent differentiation between the bones of L. 
chalconotus and L. carunculatus, Worthy [10] suggested that these 
taxa should be treated as synonyms, with regional morphological 
variation interpreted as clinal. Worthy [10] showed that the size 
range of L. carunculatus overlapped with L. chalconotus from Otago. 
Worthy [10] further stated that the only difference between the 
two taxa was the dimorphic plumage in populations of L. 
chalconotus, with the pied morph very similar to L. carunculatus, 
although, he did not consider differences in facial colour. On the 
basis of these findings and initial genetic conclusions, Schuckard 
[2 1] considered it unlikely that L. chalconotus, L. carunculatus and L. 
onslowi would maintain their specific rank in future taxonomic 
revisions. However, subsequent genetic analyses of the Phalacro- 
coracidae by Kennedy and Spencer (unpublished data) indicate 
strongly that L. carunculatus represents a distinct species that is not 
closely related to either L. chalconotus or L. onslowi. 

In light of distributional, morphological and behavioural data, a 
detailed molecular reappraisal of L. chalconotus systematics seemed 
overdue. We use both morphological and genetic data to reassess 
the systematic status of Leucocarbo shags from southern New 
Zealand. 

Materials and Methods 

Ethics statement 

Animal ethics permits from the University of Otago's Animal 
Ethics Committee were not required as samples were only taken 
from deceased birds. Specimens were not collected from 
conservation gazetted or private land, and were provided to 
NJR and MK under Department of Conservation (DoC) 
conservation management auspices. Specimens are held at the 



University of Otago's Department of Zoology under a permit to 
hold material from protected wildlife for genetic analysis (DoC 
OT-25557-DOA). Permission to sample museum specimens was 
obtained by NJR and MK from individual institutions (see 
Table S 1 for location and accession details of specimens utilised 
in this study). 

Source of specimens 

Dead, beach-wrecked L. chalconotus specimens and additional 
dead birds (including historical skins) retrieved from breeding 
colonies or roosting sites across the geographic range of L. 
chalconotus (n = 43, comprising 14 dark-bronze, 18 pied and 11 
unrecorded morphotype) were obtained from a variety of sources 
including DoC and museum collections (Fig. 1; Table SI). L. 
onslowi samples, comprising tissue and Holocene fossil bones 
(n= 10), were also included to help clarify the status of L. onslowi 
(Fig. 1; Table SI). Samples from sub-Antarctic taxa L. colensoi 
(n= 1) and L. campbelli (n = 1) were included as out-groups (Table 
SI) (see [22]). 

Modern DNA extraction, PCR amplification and DNA 
sequencing 

Whole genomic DNA was extracted from 2 mm 2 tissue samples 
following a modified Chelex protocol with 5% Chelex solution, 
2 |J.L Proteinase K (20 mg/ mL) and overnight incubation at 56°C 
[23]. For recent bone samples, DNA was extracted from up to 
200 mg of bone powder using the Qiagen DNeasy Tissue Kit 
following the manufactures instructions with an overnight 
incubation at 56°C. DNA from recent bone samples was amplified 
following the methodology for ancient DNA below. 

PCR primers were designed to amplify control region one (CR1) 
of the Phalacrocoracidae mitochondrial DNA (mtDNA) genome. 
In phalacrocoracids and core pelecaniforms there is a duplication 
of the 3' end of cytochrome b {cyt b) to the 3' end of the control 
region (CR) [24—25]. In L. chalconotus there is very low (<20%) 
sequence similarity between CR1 and CR2 compared with other 
core pelecaniforms (see [25]). CR1 was amplified using one of 
several alternate primer pairs: AV16531FND6 (5' ACCACCAR- 
CATHCCCCCYAAATA 3'; Gillian Gibb pers. comm.)/New- 
CytB Rl (5' ACAGTTTGATGAAATACCTAGTGGG 3'; this 
study) (-1200 bp including primers); AV16531FND6/NewCytB 
R2 (5' GATTTTGTCACAGTTTGATGAAATACC 3'; this 
study) (-1200 bp including primers); and AV16758FtGlu (5' 
TRTGGCYTGAAAARCCRTCGTTG 3'; Gillian Gibb pers. 
comm.)/NewCytB R2 (—1000 bp including primers). For one 
sample (Museum of New Zealand Te Papa Tongarewa (NMNZ) 
OR. 17326) only a partial CR1 fragment could be amplified using 
the primer pair AV16531FND6/H10 (5' GTGAGGTGGAC- 
GATCAATAAAT 3'; [26]). Each PCR reaction (10 uL) consisted 
of: 5 |xL of MyFi DNA Polymerase mix (Bioline), 0.5 uL each 
primer (10 |jM), and 1 |xL DNA. PCR thermocycling conditions 
consisted of 94°C 3 min, 40 cycles 94°C 30 s, 50°C 45 s, 72°C 
2 min 30 s, followed by a final extension of 72°C 4 min. PCR 
products were run on a 1 % 1 x TAE agarose gel. All extractions 
and PCRs included negative controls. PCR products were purified 
using an Omega Ultra-Sep PCR clean-up kit, then sequenced bi- 
directionally using Big Dye Terminator technology and an ABI 
3730x1. 

Ancient DNA analysis 

All ancient DNA (aDNA) extractions and PCR set up was 
carried out at the University of Otago in purpose built aDNA 
(Holocene fossil samples) and historical DNA (historical specimens) 
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laboratories physically isolated from other molecular laboratories 
[27]. Strict aDNA procedures were followed to minimise 
contamination of samples with exogenous DNA [28] including 
the use of negative extraction and PCR controls. 

DNA was extracted from up to 250 mg of bone powder 
sampled from Holocene fossil L. onslowi bones following [29] . DNA 
was extracted from L. chakonotus toe pads using the Qiagen 
DNeasy Blood and Tissue Kit following the manufactures 
instructions with an overnight incubation at 56°C. Two overlap- 
ping fragments of the mtDNA CR1 were amplified using the 
primers BESCR1 Fl (5 ' -GCCACATGATAC ATTACATG-3 ')/ 
Rl (5'-CRCTTATACATAAACTCCTAG -3') (197 bp including 
primers) and BESCR1 F2 (5 ' -C ATGTAC ARAC C C A- 
TYCCTCCC-3')/R2 (5'-GTATCCGGTTTCTGAAGTAC- 
CAG-3') (210 bp including primers). Each PCR reaction (20 |iL) 
consisted of: 1 M Betaine (Sigma), 4 mM MgCl 2 (Life Technol- 
ogies), 1 x Gold Buffer II (Life Technologies), 2.5 mM dNTPs 
(Bioline), 250 nM each primer, 1.25 U of AmpliTaq Gold DNA 
Polymerase (Life Technologies), and 2 |J.L DNA. Unsuccessful 
PCR's were repeated with 2 U AmpliTaq Gold DNA Polymerase 
and 4uL DNA or 2 uL 1:10 DNA. PCR thermocycling 
conditions consisted of 94°C 9 min, 60 cycles 94°C 30 s, 50°C 
45 s, 72°C 1 min, followed by a final extension of 72°C 10 min. 
PCR amplification and all downstream procedures were carried 
out in a modern genetics laboratory. 

PCR products were run on a 2% 1 x TAE agarose gel. PCR 
products were purified (using ExoSap (1.5 U Exol, I U SAP; GE 
Healthcare) by incubation at 37°C for 30 min and 80°C for 
15 min), and sequenced as above, except that bi-directional 
sequencing was conducted from independent PCR products. 
When an inconsistency between sequences from an individual was 
observed due to DNA damage (C-T and G-A transitions), 
additional PCRs and bi-directional sequencing were conducted, 
and a majority rule consensus was applied to the independent 
replicates [30]. 

Phylogenetic analysis 

Contiguous sequences were constructed using Sequencher 
(Gencodes) and aligned in MEGA 4.0 [31]. The alignment was 
trimmed to include 1040 bp of CR1 sequence data (i.e. excluding 
the ND6 and tRNA Glu sequence and a short region at the 3 ' end of 
CR1). ModelTest was used to determine the most appropriate 
model of nucleotide substitution under the Akaike Information 
Criterion. Two datasets were used for phylogenetic analysis: (1) 
1040 bp CR1 (modern samples) and (2) 248 bp CR1 (modern, 
recent, historic and ancient) common to all samples (for both of 
these datasets ModelTest determined the most appropriate model 
of nucleotide substitution was HKY + I). Bayesian analysis was 
performed using BEAST vl.7.4 [32] and the Yule birth-death 
coalescent tree prior. Three independent runs were conducted, 
consisting of 30 million MCMC generations, sampling tree 
parameters every 1000 generations, with a burn in of 25%. 
Phylogenetic convergence was assessed in Tracer, and results 
analysed in FigTree. DNA sequences are accessioned in GenBank 
(KJ189963-KJ190017). As well as the Bayesian posterior proba- 
bilities (PP), the level of support for the tree topology was evaluated 
with 1000 bootstrap replicates. Maximum likelihood bootstrap 
analysis was performed using PhyML [33]. Maximum parsimony 
(with equal weights and 10 random addition sequence replicates) 
and neighbour-joining bootstrap analyses were performed using 
PAUP [34]. 



Morphometric analysis 

Total length measurements (defined by [35]), of coracoids, 
humeri, ulnae, carpometacarpi, femora, tibiotarsi and tarsometa- 
tarsi were performed using vernier callipers, to the nearest 
0.1 mm. The species measured included: L. colensoi (n= 15) and 
L. chakonotus from the Otago (n = 32) and Foveaux Strait (n = 8) 
populations. These bones were sourced from recent and Holocene 
fossil skeletons housed in New Zealand museum collections 
(Table S2). Box plots of element length for each species/ 
geographic region were constructed using the program R [36]. 
Differences in physical size between L. chakonotus from Otago and 
Foveaux Strait, and L. chakonotus and L. onslowi, were assessed for 
each measured element individually, using the Mann- Whitney U 
(MWU) test, implemented in R using the wilcox.test function 
(Table 1). The MWU test is a non-parametric t-test that compares 
median values of element length distributions to determine if birds 
from one population are larger than those from another 
population. Morphological differentiation was assessed using 
Principal Component Analysis (PC A) on pooled element lengths 
(from specimens with no missing data) for L. chakonotus from Otago 
and Foveaux Strait, and L. onslowi, using the prcomp function in R. 

Results 

DNA sequence alignments for the modern specimens (15 in- 
group individuals) yielded up to 1040 bp of mtDNA CR1. Of the 
1040 bp alignment 999 bp were constant, and of the 41 variable 
sites 27 were parsimony informative. For the 248 bp mtDNA CR1 
fragment (53 in-group individuals) 218 bp were constant, and of 
the 30 variable sites 25 were parsimony informative. Our 
phylogenetic analyses of these two datasets strongly supported 
there being two clades of Leucocarbo shag in southern New Zealand 
(Figs. 2, 3, SI). The first clade comprises dark-bronze and pied 
morphotypes from Foveaux Strait along with a single dark-bronze- 
morph beach-wrecked specimen collected from the western shore 
of Boulder Beach, Otago Peninsula (collected on 22 June 2011 by 
G. Loh; currendy held at the University of Otago's Department of 
Zoology and will be deposited in the NMNZ collections). The 
second clade, which in most of the analyses is sister to the 
Chatham Islands endemic L. onslowi, comprises dark-bronze and 
pied morphotypes sampled from Otago populations. With the 
exception of maximum likelihood and neighbour joining for the 
248 bp dataset, our phylogenetic analyses produced support for 
the sister group relationship between the L. chakonotus Otago 
lineage and L. onslowi. In contrast, maximum likelihood and 
neighbour joining bootstrapping of the 248 bp dataset produced 
weak bootstrap support (54% and 58% respectively) for the sister 
group relationship between the Otago and Foveaux Strait clades 
(while still supporting the monophyly of each of these groups) 
(Fig. S2). 

For the 248 bp mtDNA CR1 fragment, the mean percentage 
divergence between the L. chakonotus Foveaux Strait and Otago 
lineages is 2.6%, whereas the Otago lineage is only 1.8% divergent 
from L. onslowi (compared to 3.5% for the contrast between the 
Foveaux Strait lineage and L. onslowi). The percentage divergence 
decreases when the longer 1040 bp of mtDNA CR1 is considered 
(L. chakonotus Foveaux Strait versus Otago 1.2%, Foveaux Strait 
versus L. onslowi 1.8%, Otago versus L. onslowi 1.0%). This 
decrease in mean percentage divergence is because a greater 
proportion of variable nucleotide positions are contained within 
the 248 bp CR1 fragment (-12% compared with ~4% in the 
1040 bp CR fragment). 

Morphometric analysis shows that L. chakonotus shags from 
Otago are larger on average than those from the corresponding 
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Table 1. Mann-Whitney U test statistics (U) and p values (P) for assessing size differentiation between Leucocarbo chalconotus from 
Otago and Foveaux Strait, and L chalconotus and L. onslowi. 
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0.00063 
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0.00049 
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<10~ 9 


Tarsometatarsus 


272.00 


0.00010 


39.00 


<10~ 9 
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Foveaux Strait populations for all elements measured (MWU P< 
0.01; Fig. 4, Table 1). The analysis also showed that, again for all 
elements, L. onslowi is smaller on average than L. chalconotus (MWU 
P<0.01; Fig. 4, Table 1). A PCA analysis of pooled element 
lengths (from specimens with no missing data) revealed substantial 
morphological differentiation between the three regional popula- 
tion groupings, with just a few immature Otago and Foveaux 
Strait birds showing morphological overlap with specimens from 
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Figure 2. Phytogeny of Leucocarbo chalconotus and L. onslowi 
shags based on 1040 bp mtDNA CRT. Branch lengths are 
proportional to the number of substitutions. The maximum clade 
credibility tree was generated in BEAST v1.7.4 using the HKY + I model 
of nucleotide substitution. The phylogeny was rooted using the 
Auckland Island Shag (L colensoi) and Campbell Island Shag (L 
campbelli). For clarity, only posterior probability (PP) and bootstrap 
support (bold: maximum likelihood; italics: maximum parsimony; 
Roman: neighbour joining) values for major clades are shown (>0.60 
PP). Nodes with less than 0.60 PP have been collapsed. Dark-bronze 
(black circle) or pied (white circle) plumage morphotype, if known, has 
been indicated on the phylogeny. 
doi:1 0.1 371 /journal.pone.0090769.g002 



the Chatham Islands (Fig. 5). Even with these immature birds 
included, the morphometric differences observed between L. 
chalconotus from Otago and Foveaux Strait, and L. chalconotus and L. 
onslowi, are highly significant. 

Discussion 

Our phylogenetic analyses show that Leucocarbo shags in 
southern New Zealand comprise two well-supported clades (Otago 
and Foveaux Strait), each containing both pied and dark-bronze 
morphs. Surprisingly, however, the combined monophyly of these 
L. chalconotus populations is not well supported, with the Otago 
lineage typically found to be sister to the Chatham Island endemic, 
L. onslowi, albeit with weak statistical support (a small subset of our 
analyses do weakly support the monophyly of the two L. chalconotus 
populations). Morphometric analysis indicates that Leucocarbo shags 
from Otago are larger on average than those from Foveaux Strait, 
corroborating Lalas' [4] observations of substantial size differences 
between individuals from these regions. 

We interpret the specimens from the Southland coast (Oreti 
Beach and Invercargill) as beach-wrecked individuals from the 
Foveaux Strait population (Figs. 1-3, SI). Given the prevailing 
southerly winds in the region, it is entirely plausible that dead 
individuals from Foveaux Strait will regularly wash up on southern 
South Island beaches. We also interpret the beach-wrecked 
specimen from Boulder Beach, Otago Peninsula, as likely of 
southern origin. This hypothesis is supported by the specimen 
having scattered papillae, consistent with the Foveaux Strait 
morphotype [4] (Fig. 1). Juvenile and adult Leucocarbo shags are 
known to travel long distances [6], [13], which, in combination 
with the direction of the Southland Current, suggests that finding a 
beach-wreck from the Foveaux Strait population on the southern 
side of Otago Peninsula is not altogether unexpected. Although we 
did not sample specimens from every location that the species 
occur, from external morphology we would also expect extant 
individuals from Nugget Point, Kinakina Island and The Sisters in 
the Catlins (plumage proportions reported by [5] are 20-30% pied 
for these colonies) to represent the Otago lineage (Fig. 1). The 
specimen from Cooks Head, Milton, just north of Nugget Point, 
falls within the Otago clade, which is consistent with this 
suggestion. 

Our phylogenetic results, together with consistent differences in 
relative proportions of plumage morphs and facial carunculation, 
and concordant differentiation in body size and breeding time [4— 
5], suggest three taxonomic alternatives requiring further inves- 
tigation that are all consistent with the observed phylogeographic 
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Figure 3. Phytogeny of Leucocarbo chalconotus and L. onslowi shags based on 248 bp mtDNA CR1. Branch lengths are proportional to the 
number of substitutions. The maximum clade credibility tree was generated in BEAST v1.7.4 using the HKY + I model of nucleotide substitution. The 
phylogeny was rooted using the Auckland Island Shag (L colensoi) and Campbell Island Shag (L campbelli). For clarity, only posterior probability (PP) 
and bootstrap support (bold: maximum likelihood; Italics: maximum parsimony; Roman: neighbour joining) values for major clades are shown (>0.60 
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PP). Nodes with less than 0.60 PP have been collapsed. Weakly supported alternative branching patterns are not shown. Dark-bronze (black circle) or 
pied (white circle) plumage morphotype, if known, has been indicated on the phylogeny. ': 0.88, 55, 53; 2 : 1 .00, 75, 78, 82; 3 : 0.88, 62, 55, 65; 4 : 0.98, 
68, 59, 73. 
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pattern of (Foveaux (Otago, Chatham Islands)) (Figs. 2, 3, SI). 
Resulting taxonomic classifications depend on differing interpre- 
tations of taxonomic philosophy. The first hypothesis is that the 
Otago and Foveaux Strait lineages of L. chalconotus represent 
phylogenetically and diagnosably distinct taxa (at the species level), 
rather than the prevailing view that there is a single variable taxon 
[1—2], [10]. This hypothesis also recognises the separate specific 
status of the rare Chatham Island endemic, L. onslowi. L. onslowi 
individuals are morphologically distinct as they are only of the pied 
morph and smaller on average than L. chalconotus (Figs. 4, 5), and 
their ranges do not overlap [1], [14]. Worthy [20] found 14 



Cora 



Ulna 



Fem 



Hum 



Cmc 



Tt 



Tmt 



Chathams Otago Foveaux 



Chathams Otago Foveaux 

Figure 4. Box plots showing the distribution of Leucocarbo 
chalconotus and L. onslowi element lengths (mm). From the 
median (thick black line) outwards, are the 75th and 25th percentiles 
(box), the 98th and 2nd percentile (whiskers) and outliers (hollow dots). 
L chalconotus from Otago are larger on average than Foveaux Strait 
birds (P<0.01; Table 1), and L. chalconotus are larger than L. onslowi 
(Table 1). 

doi:1 0.1 371 /journal.pone.0090769.g004 



osteological character state differences between L. chalconotus and 
L. onslowi. L. onslowi individuals also differ substantially in facial 
colouring from L. chalconotus. Breeding L. onslowi individuals 
consistently have pronounced bright orange caruncles compared 
to small bright orange caruncles (Otago) and scattered dark to dull 
orange papillae (Foveaux Strait and Otago) in L. chalconotus [13], 
Though more variable, but following a similar pattern to the 
carunculations, the gular throat patch is bright red in breeding L. 
onslowi individuals, while it ranges from dull to bright orange-red in 
L. chalconotus [13]. Determining whether the Otago and Foveaux 
Strait lineages of L. chalconotus represent diagnosably distinct taxa 
will require (1) osteological analyses of recent skeletons (identified 
to morph, age and sex) and Holocene fossil bones to test for the 
presence of discrete morphological characters; (2) additional 
morphological analysis of the extent of carunculation to determine 
whether this can be used as a discrete morphological character; (3) 
genetic analysis of specimens from un-sampled areas (e.g. 
Kinakina Island, The Sisters in the Catlins); and (4) ancient 
DNA analyses of type specimens, Holocene fossil and archaeo- 
logical material. 

The second plausible hypothesis is that L. chalconotus may 
represent a paraphyletic species with strong phylogeographic 
structure, a consequence of the regional philopatry that is well- 
documented in this group of birds. Such a finding would not be 
unexpected, as Joseph and Omland [37], for instance, showed that 
44% of Australo-Papuan terrestrial avian taxa that are good 
biological species are paraphyletic for mtDNA markers. Applying 
a biological or phylogenetic species concept to taxa with allopatric 
populations can be problematic, especially when they exhibit site- 
specific or regional philopatry. Future research will focus on 
population-genetic analysis of the Otago and Foveaux Strait clades 
to determine the level of gene flow, population divergence and 
dynamics between each region [38-39]. There are only two single 
nucleotide polymorphisms (out of over 5000 bp of nuclear DNA) 
that separate L. onslowi and the L. chalconotus Otago lineage 
(Kennedy and Spencer unpublished data). As a consequence, 



o 
O 



■ Otago 

■ Foveaux Strait 
Chatham Is. 



Principal Component 1 

Figure 5. PCA cluster analysis of Leucocarbo chalconotus from 
southern New Zealand, and L. onslowi. The analysis was conducted 
on pooled element lengths (femur, tibiotarsus, tarsometatarsus, 
coracoid, humerus, ulna, and carpometacarpus). Immature birds are 
indicated with an T on the figure. 
doi:10.1371/journal.pone.0090769.g005 
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genetic discrimination of the Otago and Foveaux lineages based 
on nuclear sequences is unlikely. 

As a consequence of this second hypothesis, a third hypothesis 
could be advocated that L. chakonotus and L. onslowi represent a 
single, monophyletic taxon with marked phylogeographic struc- 
ture and possible sub-specific status for the Foveaux, Otago and 
Chatham Island clades. Nevertheless, this hypothesis is unlikely as 
L. onslowi is morphologically (Figs. 1, 4 and 5) [1], [14], [20] and 
genetically (Figs. 2, 3, SI) distinct from L. chakonotus. 

The species-rich radiation of New Zealand blue eyed shags 
(including L. ranfurlyi, L. campbelli, L. colensoi, L. onslowi and L. 
carunculatus), and especially the diversity within L. chakonotus, 
contrasts strongly with the situation for South American/ Antarctic 
blue eyed shags which show less evidence for evolutionary 
diversification [1]. The phylogeographic partitioning detected 
within L. chakonotus in particular is marked, and such strong 
structure is rare for Phalacrocoracidae species. Within Leucocarbo, 
Calderon et al. [39] found strong phylogeographic structure and 
limited gene flow in the Patagonian Rock Shag (Leucocarbo 
magellanicus) using mtDNA ATPase 6-8 and microsatellites. The 
Rock Shag contained three mtDNA phylogeographic groups 
comprising colonies from the Atlantic, Pacific and Fuegian regions 
of Patagonia. In contrast, there was weaker phylogeographic 
structure and higher levels of gene flow in the sympatric Imperial 
Shag (Leucocarbo atriceps) [39]. Calderon et al. [39] concluded that 
Pleistocene glacial cycles, and differences in foraging ecology and 
non-breeding distribution between the Rock and Imperial Shag 
contributed to the differing levels of phylogeographic structure. 
Winney et al. [26] and Marion and Gentil [40] found only 
moderate phylogeographic structure in the Great Cormorant 
(Phalacrocorax carbo) from Europe using mtDNA CR1 data, while 
Barlow et al. [41], using a combination of mtDNA ND2 and 
microsatellite data from the European Shag (Phalacrocorax aristotelis), 
found only weak population-genetic structure, despite high levels 
of philopatry, rare dispersal events and recognised subspecies. 
Waits et al. [42], using mtDNA 12S rRNA, 16S rRNA, COIII-ND3, 
cyt b and CR, found no phylogeographic structure or restricted 
gene flow among populations or between subspecies of the 
Double-crested Cormorant (Phalacrocorax auritus) in eastern North 
America. However, Mercer et al. [43] used CR data and found 
there was phylogeographic structure between the Alaskan Double- 
crested Cormorant and those from eastern North America. In 
another coastal bird taxon, the Dusky Seaside Sparrow (Ammo- 
dramus maritimus), Avise and Nelson [44] detected substantial 
phylogeographic structuring of populations from eastern and 
southern North America, delineated by the Florida peninsula. It 
should be noted that the genetic data used in our study are from 
rapidly evolving mtDNA CR1, rather than the more slowly 
evolving 'barcoding' COI marker often used in bird systematics 
(e.g. [45] cf. [46]). However, there are several diagnosably distinct 
avian taxa than cannot be distinguished using slowly evolving 
mtDNA markers but can be identified using the faster evolving CR 
(e.g. cyt b versus CR in Cyanoramphus parakeets [47-48]). 

Conclusions 

In this study we have shown that the Stewart Island Shag (L. 
chakonotus) groups into two geographically separate (Otago and 
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